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Abstract 

In this paper an analytical and numerical study of anharmonic vibrations of 
monatomic chain and graphene in transverse (perpendicular) with respect to the 
chain/plane direction is presented. Due to the lack of odd anharmonicities and 
presence of hard quartic anharmonicity for displacements in this direction, there 
may exist localized anharmonic transverse modes with the frequencies above the 
spectrum of the corresponding phonons. Although these frequencies are in resonance 
with longitudinal (chain) or in-plane (graphene) phonons, the modes can decay only 
due to a weak anharmonic process. Therefore the lifetime of these vibrations may be 
very long. E.g. in the chain, according to our theoretical and numerical calculations 
it may exceed 10^^ periods. We call these vibrations as transverse intrinsic localized 
modes. 
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1 Introduction 


It is a well-known fact that point defects in crystals may canse an appear¬ 
ance of spatially localized vibrations, called local modes (see, e.g. m)- The 
freqnency of such modes lies outside the phonon spectrum, which prevents 
the spreading of the vibrations to the bulk. Besides, resonant or pseudolocal 
modes with the frequencies inside the phonon spectrum may also exist. They 
appear when the frequency of the leading vibration associated with the de¬ 
fect gets to a region of small density of states (DOS) of phonons. Unlike local 
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modes, which are stable in harmonic approximation, these modes can live only 
a finite time due to the emission of resonant phonons. However, because of 
the small DOS of resonant phonons the interaction of the latter with pseu¬ 
dolocal mode is weak, due to what the lifetime of the mode becomes long. A 
well-known case of long-living pseudolocal modes in 3D lattices is the one by 
the defect of large mass [T]. 

Spatially localized modes may exist also in perfect nonlinear lattices. Such 
modes in chains with cubic and quartic anharmonicity were first described 
by A. M. Kosevich and A. S. Kovalev [2]. They have found that for suffi¬ 
ciently strong quartic anharmonicity there exist vibrational states localized 
in the space and periodic in time; the frequency of these vibrations exceeds 
the maximum frequency of phonons. The authors [2] restricted their consid¬ 
eration with large-size vibrations and the frequency close to the top of the 
phonon spectrum. Small size localized vibrations in anharmonic lattices were 
introduced by Dolgov |3] and then by A. J. Sievers and S. Takeno [1] and were 
called “intrinsic localized modes” (ILMs). 

In the numerical studies of ILMs different two-body potentials (Lennard-Jones, 
Born-Mayer-Coulomb, Toda, and Morse potentials and their combinations) 
have been used (see, e.g. [51I6H7] ). All these potentials have strong odd anhar- 
monicities and show a strong softening with the increasing of the vibrational 
amplitude. The ILMs found in these simulations always drop down from op¬ 
tical band(s) into the phonon gap, if there is any. In this connection see Ref. 
[S], where gap ILMs in Nal were calculated without taking into account long- 
range interactions and Ref. [9], where the calculations of ILMs in Nal were 
made taking these interactions into account. 

Usually in crystal lattices odd anharmonicities are strong and pair potentials 
show a strong softening with the increase of vibrational amplitude. Therefore, 
the dropping down of the frequency of ILMs from optical bands is quite com¬ 
mon in 3D lattices. Still, as it was shown in Refs. [T01I11II12I1T3] . in some crystals 
odd anharmonicities are reduced due to multiparticle or covalent interactions. 
Examples have been given with germanium [10], diamond m, metallic Ni, 
Nb [12], iron [TB] and copper m]; in all these crystals ILMs with frequencies 
above the top of the phonon band were found in numerical simulations. 

However, there are systems in which the odd anharmonicities disappear due to 
symmetry arguments; at the same time, the quartic anharmonicity is non-zero 
and it is hard (positive). The examples of such systems are given by linear 
atomic chains and planar atomic structures (e.g. graphene); in these systems 
odd anharmonicities disappear for vibrations in the transverse (out-of-chain 
and out-of-plane) direction. Therefore, one can expect that in these systems 
there can exist transverse anharmonic local modes with the frequencies above 
the maximum frequency of the corresponding phonons. The latter frequen- 
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cies are usually smaller than the maximum frequency of longitudinal/in-plane 
phonons. These modes fall in resonance with these phonons and can decay. 
However, unlike pseudolocal modes and like local modes the interaction caus¬ 
ing the decay of these modes is anharmonic. Therefore we call them as trans¬ 
verse ILMs (TILMs). The anharmonic interaction of TILMs with small vi¬ 
brational amplitudes of atoms is very weak. Hence, one can expect that the 
lifetime of such TILM may be rather long. Below we will consider TILMs in 
a monatomic chain and in graphene, both analytically and numerically. Our 
considerations conhrm the aforesaid expectation. 


2 Anharmonic chain 


First we consider the anharmonic monatomic chain and examine the vibrations 
of its atoms in transverse {y) directions. We suppose that the potential energy 
U of the chain is given by the sum of pair potentials V{Rn,n'), where Rn,n' is the 
distance between the atoms n and n'. The latter potentials can be expanded 
into the series of atomic displacements. Denoting d = {n — n') a, x = Xn — Xn', 
y = Vn — Un'i where a is the atom spacing, and yn are the longitudinal and 
transverse displacements of the atom number n from its equilibrium position, 
we get R = Rn,n' = \J{d + xY + y'^. As y appear as y'^, any power expansion 
of R will have only even powers of y. The same holds for U. This means that U 
indeed has no odd anharmonic terms. This is a consequence of the symmetry 
of the chain with respect to the change of the sign of y. 

Note one more property of the chain; the term in the expansion of R, quadratic 
with respect to y, has the same numerical factor as the term linear with respect 
to X. The same holds also for an arbitrary power of R. In the equilibrium state 
all linear terms with respect to the coordinates x in the potential energy U are 
cancelled. Therefore, all quadratic terms with respect to y are also cancelled, 
i.e the frequencies of transverse vibrations in the pair potential approximation 
tend to zero. As a result long-range fluctuations can be created with little 
energy cost and since they increase the entropy they are favored. This leads 
to the instability of the chain with respect to small transverse distortions (see 
in this connection the Mermin-Wagner theorem mi)- 

To get the chain stable one needs to stretch it [15]. In this case the atom 
spacing a is replaced by a -|- s, where s is stretching. Then the terms in U 
linear with respect to x and are not cancelled any more. Therefore, the 
elastic springs for transverse vibrations are also nonzero and positive. This 
results in the appearance of transverse phonons with hnite, although small 
for small stretching maximum frequency ujtm ■ As these phonons do not have 
any cubic anharmonicity, but have nonzero positive quartic anharmonicity, 
low-frequency TILMs with the frequency above the spectrum of transverse 
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phonons should exist here. 

Let us consider the TILM in a monatomic chain with the Morse pair potential 

= ( 1 ) 


Here D is the energy of dissociation, a is the parameter. We are using dimen¬ 
sionless coordinates with the units corresponding to a = 1 and the value a = 4 
of the Morse pair potentials of atoms in monatomic metals. For this poten¬ 
tial only nearest-neighbor interactions are essential and only these potentials 
will be taken into account here. We also take for the mass units the mass of 
the atoms of the chain (M = 1). The dissociation energy is chosen so that 
the unit frequency will correspond to the maximum frequency of longitudinal 
phonons. In this case the potential energy of the stretched lattice is the sum 
of the following pair potentials (up to a constant term): 

1 / = (l - y 128 - a: (l - 16, ( 2 ) 


where r = -|- 1 -f s)^ -|- The last term in Eq. (E]) accounts for the effect 

of the stretching force of the chain in x direction - it changes the equilibrium 
distance of the atoms from r = ltor = l-|-s. Let us expand the potential 
into the series of x and y and take into account up to the second-order terms 
with respect to x and forth-order terms with respect to y. We get (up to a 
constant) 


T// ^ ^ 2 I ^2 2 I ^3 2 I ^4 4 

V{x,y)^—x +—sy +—xy + ^2/ ) 


( 3 ) 


where Ui are dependent on stretching s parameters. In the small s limit i/j ^ 1. 
If s = 0.05 then z/i 0.522, 1/2 ~ 0.707, z /3 ^ 0.463, z /4 ^ 0.441. 

The pair potential of the longitudinal vibrations alone is given by the hrst 
term in the right-hand side of Eq. ([3]). Vibrational frequencies of corresponding 
phonons equal [T] 




cos{k)) / 2 . 


( 4 ) 


The maximum frequency of longitudinal phonons corresponds to k = —it and 
equals uim = \/fq- T^e transverse vibrations alone are described by the pair 
potential 

Vi0,y) = ^y‘^ + ^y^. (5) 
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In harmonic approximation (z /4 = 0) the frequencies of corresponding phonons 
are given by Eq. (|1]) with 1 / 2 S instead on z/i. The positive quartic anharmonic- 
ity in Eq. (jS]) leads to appearance of the anharmonic modes [2] (called here 
as TILMs) with the frequencies ujq = above the maximum fre¬ 

quency of transverse phonons Utm = y/^ 2 S and with the displacements 

Unit) ~ cosh"^(£n) cos(a;ot) . ( 6 ) 

Here Aq is the amplitude of the central atom, 

(7) 


is the reversed size of the TILM (we use the discrete analog of the derived 
in Ref. |2] equation (47) for the difference in the displacements of two neigh¬ 
bouring atoms x(^) )• These modes interact with longitudinal phonons and, 
therefore, they decay. To describe this decay we consider the longitudinal vi¬ 
brations of atoms in the presence of the TILM. Taking into account Eq. (jH]) 
we replace y by y{t) in Eq. (|3]) and get the following pair potential for this 
motion: 

V(x, t) ^ + ^x|/^( 0 ) -f- ^xy‘^{fi) cos( 2 cnot). ( 8 ) 

8 lo lo 


The hrst term in the right-hand side of this equation gives the potential energy 
of the longitudinal vibrations alone in harmonic approximation. The second 
and the third terms describe the anharmonic interaction of these and trans¬ 
verse vibrations. At that the second term stands for a small local compression, 
while the third term describes the force with the frequency 2ujq periodically 
changing in time. For all longitudinal phonons, except those with the resonant 
frequency Uk = 2uo, this force causes forced vibrations of atoms with the fre¬ 
quency 2uo- The resonant term causes the increase of the energy of phonons 
in time. From energy conservation law it follows that this energy comes from 
the TILM, i.e. the TILM decays. 


3 Decay rate of transverse ILM in chain 


To hnd the rate of decay we are considering the equation of motion of the 
longitudinal phonon with the coordinate 

Xk = Xn sin(/cn), (9) 
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where S> 1 is the number of atoms in the chain, k = 27rm/iV is the wave 
number of the phonon, m = 0, ±1, ib2, ... (we use periodic boundary condi¬ 
tions), Xn is the displacement of the atom n from its equilibrium position. In 
the case of pair potential (|H]) this equation reads 

Xk = -ojlxk + sin(fcn)(^(2/o,n+i - yo,nf - ( 2 / 0 ,n - yo,n-if) cos(2a;ot), (10) 


where Uk is the frequency of the longitudinal phonon k given by Eq. (|1]), 

= 2/n(0) is the initial amplitude of the transverse ILM on atom n given 
by Eq. (E]). The hrst term in the right-hand side of Eq. fllOp stands for the 
harmonic force stemming from the hrst term in the right-hand side of Eq. 
(18|) and the second term accounts for the anharmonic force stemming from 
the third term in the right-hand side of this equation (a small compression of 
the chain stemming from the second term in the right-hand side of Eq. ([8]) is 
neglected). The transverse amplitude |2/o,n| slowly changes with n. Therefore, 
2 / 0 ,n+i 2 / 0 ,n — 2 j/q ,1 and 

( 2 / 0 ,n+l - 2/0,n)^ - ( 2 / 0 ,n “ 2/0,n-l)^ ~ 4 dy^/du . (11) 


Inserting Eq. fllip into Eq. m and taking into account Eq. (E]), we get 

Xk +ulxk = CkN~^^'^ cos(2a;o t), (12) 


where 


Ck = ^Al Y sin(/cn 



cosm(£n) ^ 


= —eusA-l Y sm(fcn) 


n=0 


sinh(£n) 

cosh^(£n) 


.(13) 


By using the Green’s function of harmonic oscillator Gk{t) = ujk^ sm{ujkt), the 
solution of Eq. (IT^ can be presented in the following form: 

Xk{t) = Xo^k cos{ukt) + x'^{t), (14) 


where 

C r 

x'f,{t) = — / sin {ujk{t— t'))cos{2ujot') dt'. (15) 

^/Nuk I 

The hrst term in Eq. flTTp (oc XQ^k) stands for free oscillations, while the second 
term (a;).(t)) describes the excitation of vibrations by the transverse mode. 
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The energy of the excited vibrations is given by the snm of the terms x'\ 
averaged over period and snmmed over all phonons. This gives 

1 * * 

E{t) = dt'sin {uk{t — t')) cos {2uot') j dt''sm. {ujk{t — t")) cos (2a;ot^0)(16) 

^ 0 0 

Let ns take into acconnt the relations: 

sin(a;A;(t — t')) sin(a;fc(f — t”)) = (cos(ci;fc(T — t")) — cos{{uJk{2t — t’ — t")))/2 , 
cos(2ci;ot') cos(2ci;ot") = (cos(2ci;o(t' — t")) + cos(2ci;o(t' + t")))/2 . 


In the case nnder consideration freqnencies of phonons form a continuous 
spectrum. In this case, due to the summation (integration) over k the terms 
in Eq. (flbll oc cos(ci;a;(T — t”)) and oc cos((a;fc(2f — t' — t”)) essentially differ 
from zero only for hnite values of the arguments t' — t" and 2t — t' — t" (of the 
order of unity). At that the term oc cos((ci;A;(2f — t' — t"))) essentially differs 
from zero only for large T, t" ~ t when cos(2a;ot') and cos(2a;ot") average the 
contrubution of this term out. Neglecting now unessential for T, t" ~ t » 
term oc cos(2a;o(t' + 1”))) we come to conclusion that at large t only the term 
oc cos(2a;o(t^ — t")) can give a remarkable contribution to the energy E{t). 
Finally we take into account that at large time t in the factor 


cos{uk{t'—t")) cos{2u}o{t'—t”)) = {cos{{uJk—2uo){t'—t"))+cos{{uJk+2uo){t'—t")))/2 


only the term cos((ci;o — 2uo)(t' — t"))/2 is essential. This gives 
1 


t t 


= /*'/ dt" cos ({ut - 2uia)(t'-t")) 


( 17 ) 


0 0 


t t 

For small \t' —1"\ being essential here one can replace the integral / / dt'dt” = 

0 0 

it' t oo 

f dti f dr by the integral / dti J dr, where E = {t' + 1")/2, r = t' — t". As a 

0 —t' 0 —oo 

result, in the t ^ oo limit the rate of increasing of the energy of longitudinal 
phonons, followed by a subsequent decrease of the energy of the TILM under 
consideration, equals 


dE{t)/dt = {ti/AN) ^ 0^6{uk — 2 u;q), 

k 


(18) 
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Replacing here the sum by the integral and taking into account the relations 
dm/dk = N/2 'k and dk/dujk = sink one hnds 


dE/dt K. [uqvIAq! sin/co) (^), 


(19) 



OO 



4(S1) = I 


( 20 ) 


0 


The decay constant T of the TILM is determined by the equation T = E ^dE/dt 
Taking into account the relation E = uJq I]n2/n(0)/2 = 


r = (fi), 


( 21 ) 


ojq sin fco 


T decreases with decreasing of the amplitude of the central atom Aq of the 
TILM. To get an estimation we consider the TILM with the amplitude Aq = 
0.03 in the 5% stretched chain (s = 0.05). In this chain ujtm ~ 0.1880, 
e: 0.183, ujq ~ 0.1888, which gives ^ 6 and <I>^ (fl) ~ 2 ■ 10“®. As a result 

we get T ^ 5 • 10“®. This corresponds to a very long lifetime. 


4 Transverse ILM in the chain; nnmerical study 

In our numerical study of transverse ILMs we considered a 5% stretched chain 
with 40 000 atoms and with hxed ends. The small compression of the entire 
chain due to TILM was taken into account. The potential energy was described 
by the interactions of the nearest neighbors with the Morse pair potential given 
by Eq. o. We have found that TILMs indeed exist in this chain and they 
have a long lifetime (see Figs. 1-3; one period of the TILM with the frequency 
cjq ~ 0.1888 corresponds to 33 time units). E.g. in our numerical simulations 
we have found that the TILM with Ao = 0.03 and s = 0.05 has the frequency 
uq ~ 0.1883 and the reversed size parameter e 0.113. It decays very slowly: 
its amplitude diminishes less than 2 ■ 10“® for 1000 periods of vibrations. This 
corresponds to T < 2-10“^^, which indeed gives very long lifetime, in agreement 
with the conclusion of the theory. 
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Fig. 1. Time dependence of the transverse shifts (in the units of atom spacing) of 
the central atom of a transverse IPLM in 5% stretched monatomic chain; the initial 
shift of the central atom is 0.05. First 10 periods and 10 periods after 1000 periods 
of vibration are presented. 

The TILMs with the amplitudes > 0.05 decay faster; the rate of their 
decay rapidly increases with the increasing of the amplitude. This is also in 
agreement with the theory. 



Fig. 2. The same as in Fig. 1, but the initial shift of the central atom of a transverse 
ILM is 0.15. 
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Fig. 3. Time dependence of the amplitudes of the central atom of transverse ILMs 
in 5% stretched monatomic chain for different initial shifts of the central atom. 


Concerning the difference of the above-presented numerically-found values of 
the parameters £ and uq from their theoretical values, we have found that 
these differences come from the neglecting of the local compression of the 
chain in the theory. This follows from the numerical simulations of TILMs 
in the chain rigid in x direction (with x coordinates of atoms being hxed at 
their values without TILM), performed by us. We have found the TILM in 
this chain with £ 0.183 and cjq ~ 0.1888 in full agreement with the above 

theoretically found values of these parameters. This TILM does not decay, as 
it should according to the theory presented above. An approximate account of 
the compression mentioned in the theory, allowing one to explain the values 
e ^ 0.113 and ujq ~ 0.1883 found in simulations, is given in Appendix. Note 
that if to use these values of £ and ojq in Eq. (1^ we get T 2 • 10“^^. This 
also corresponds to a very long lifetime. 
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5 Graphene: out-of-plane vibrations 


The arguments presented in the beginning of Section 2 can readily be extended 
to planar atomic lattices. To be specific, we are considering here a graphene 
sheet. It consists of carbon atoms constituting the honey comb-type period¬ 
ical structure. Every carbon atom in this structure is connected with three 
neighboring atoms by chemical bonds formed by the sp 2 hybrid orbitales. 

We suppose that in equilibrium positions the graphene is situated in the (x, y) 
plane. The distances between an atom and three neighboring atoms equal 



( 22 ) 


{n = 1, 2, 3), where Xm Vn and Zn are the x-, y- and 2 ;- components of the dis¬ 
placements of three neighboring atoms with respect to the central atom, bi = 
0, 62 = \/3/2, &3 = , Cl = —1, C2 = 1/2, C3 = 1/2. In the pair potential 

approximation the potential energy of the vibrations of atoms depends on the 
distances — a (explicit form of this energy in harmonic approxima¬ 

tion, see e.g. in [I^). The expansion of these distances over the displacements 
Xn, yn and Zn depends on the powers = {hx + cy + [x^ + y'^ + z‘^)/a)^ 
(here the subscript n is omitted for simplicity). In the equilibrium position 
the oc R terms are cancelled. Therefore, the expansion of a pair potential over 
z aX X = y = Q starts with the positive quartic term oc 2 ;^, i.e. the oc z"^ terms 
are absent. Consequently, in the pair potential approximation the 2D lattice is 
unstable with respect to small out-of-plane distortions. However, in graphene 
the atomic interactions are determined by covalent forces. These forces cause 
the stiffness of the planes with respect to the transverse short-range displace¬ 
ments of atoms [16]. In accordance with the Mermin-Wagner theorem [IT] 
the plane of graphene remains unstable with respect to long-range transverse 
distortions, resulting in the appearance of ripples [T7IIT8] . A stretching of the 
graphene sheet removes the ripples. 

The elementary cell of graphene includes two atoms. Therefore, there are two 
branches of out-of-plane vibrations: acoustic and optic. One can expect the ex¬ 
istence of two types of out-of-plane TILMs in graphene - with the frequency 
above the top frequency of acoustic out-of-plane phonons (the acoustic-like 
TILM) and with the frequency above the top frequency of optic out-of-plane 
phonons (the optic-like TILM). The acoustic-like out-of-plane TILMs can de¬ 
cay due to a relatively weak anharmonic interaction with the in-plane phonons 
and also due to in general stronger harmonic resonant interaction with out- 
of-plane optic phonons. However, the optic-like out-of-plane TILMs can decay 
only due to a relatively weak anharmonic interaction with in-plane phonons. 
Therefore, the optic-type out-of-plane TILMs should have a longer lifetime. 
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This conclusion is in agreement with our numerical modelling of out-of-plane 
TILMs in weakly-stretched (1.5 %) graphene. We have used the AIREBO 
potential and taken into account the cluster 246 x 210 A (100 x 50 periods, 20 
000 atoms) with periodical boundary conditions. Initially we displaced out-of- 
plane six atoms of Ce ring according to the optical mode. We found that out- 
of-plane optic TILMs indeed exist; see Fig. 4 (where the time-dependences of 
the displacements of one of the three equivalent atoms of the central Ce ring 
are presented; the other three atoms in this ring vibrate out of phase with 
the mentioned atom) and Fig. 5 (where in-plane and out-of-plane phonons 
DOS and the spectrum of out-of-plane TILM for different initial amplitudes 
of carbon atoms are given). 




Fig. 4. Optical-type out-of-plane transverse ILMs in graphene. The out-of-plane 
shifts (independent of time) of one of three equivalent atoms of the central Ce ring 
is presented for initial amplitudes 0.4 A (upper) and 0.2 A (bottom); the other three 
atoms in this ring vibrate out of phase with the atom mentioned. 


12 






Fig. 5. In-plane and out-of-plane phonon DOS for the AIREBO potential and the 
spectrum of out-of-plane transverse ILM for different initial amplitudes of carbon 
atoms. 


Note that in strongly-stretched graphene there can exist also in-plane anhar- 
monic localized vibrations (in-plane discrete breathers) [TM20] . It appears (see 
Ref. [19]) that strong nniaxial stretching (along the zigzag or armchair direc¬ 
tion) results in the opening of the gap in the middle of the phonon spectrum 
of graphene. This makes possible the existence of in-plane soft ILMs (discrete 
breathers) with the frequency in the gap of the phonon spectrum [T9]|^ . The 
authors [ISEQI have found that these modes remain relatively stable also when, 
due to a large stretching, their frequencies get into resonance with out-of-plane 
phonons. 


Thus, the physical situation discussed in Refs. [151120] is reversed to the one 
considered here: the anharmonic in-plane modes with soft anharmonicity are 
in resonance with out-of-plane phonons, while out-of-plain TILMs considered 
here are in resonance with in-plane phonons. However, in both cases the life¬ 
time of the modes is long. 
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6 Discussion 


In this communication it was shown that due to symmetry arguments in lin¬ 
ear monatomic chains and in planar monatomic structures (e.g. graphene) 
odd anharmonicities disappear for the motion of atoms normal with respect 
to the chain/plane. However, even (quartic) anharmonicity exists and is hard 
(positive). As a result, the a n h armonic localized modes of transverse (chain) 
or out-of-plane (graphene) vibrations (called TILMs) can exist in these sys¬ 
tems with frequencies above the transverse/out-of-plane phonons. These vi¬ 
brations have finite, although large lifetime: the slowly decay due to creation 
of longitudinal phonons. For the case of a monatomic chain a theory is being 
developed which allows one to calculate the decay rate of TILMs. The theory 
predicts that the lifetime of the TILMs under consideration may be extremely 
long, longer than 10^*^ periods of vibration. The numerical modelling of these 
TILMs fully confirms this result. The normal (out-of-plane) TILMs with the 
frequency above the phonon spectrum of out-of-plane phonons were also found 
for weakly-stretched graphene. 

Here only immobile TILMs were considered. However, one can expect that 
these excitations can move in the same way as ILMs can move in Cu, Ni and 
Fe in the nearest neighbors crystallographic directions [TT1I12II13] . 

The symmetry argument presented here holds generally for linear atomic 
chains and planar monolayer atomic structures. Therefore, one can expect 
the existence of TILMs in any linear atomic chain and any planar monolayer 
atomic structure. Besides, one can expect that an analogous situation may 
exist in 3D lattices with dominating chain and layer structural elements, e.g. 
in graphite and in mica. In such lattices odd anharmonicities exist, but they 
may be essentially reduced for out-of-chain/plane vibrations. 

In this connection we would like to point to investigations of M. Russell et 
al. of muscovite mica irradiated by high energy particles (see, e.g. [21]). The 
irradiation produces in mica very long black tracks parallel to the crystallo¬ 
graphic plains and localized between the layers. These tracks are explained in 
Ref. |2I] by creation in the recoil process of discrete breathers which propagate 
on long distance along the crystallographic directions. The layered structure 
on atomic scale plays a crucial role in the observed phenomena. Although it is 
not yet clear whether these long-living vibrational excitations are transverse 
or longitudinal, the necessity of the layered structure on atomic scale to ob¬ 
serve the phenomenon gives support to an assumption that these excitation 
should contain remarkable transverse component. 

Analogous TILMs may be expected to exist also for surface vibrations in 
crystals with planar atomic surfaces. The change of the potential energy of 
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a surface atom for the shift in outer direction also does not have any odd 
terms with respect to the shift. This leads to the decrease of the entire odd- 
anharmonicity effects for out-of-plane surface vibrations, which may result in 
the appearance of out-of-plane TILMs with frequencies above the maximum 
frequency of transverse surface phonons and with a long lifetime. 

It is known that next-to-next interactions in a chain can cause its zigzag- 
type secondary structure [22]. In such chains, e.g. in polyethylene, the mobile 
“solitons of tension” may exist [23] (involving the longitudinal-type motion of 
atoms). Besides, in such chains one should also expect the existence of trans¬ 
verse ILMs. Here two types of transverse vibrations should be distinguished: 
parallel and perpendicular to the plain of zigzag. The odd anharmonic terms 
for the transverse displacements perpendicular to the plain mentioned are 
absent; however, for transverse displacements parallel to this plain these are 
present. If the zigzag angle is large, then the TILMs perpendicular to the 
plain of zigzag should exist. However, if this angle is small, then the odd an- 
harmonicity for the displacements parallel to the plain of zigzag is reduced. In 
such chains the existence of both TILMs is expected. It is worthwhile to note 
that intrinsic localized modes (discrete breathers) are often used to explain 
the targeted energy transfer in such basic for living organism organic chains 
as DNA and proteins [2B] . 

Acknowledgements. The research was supported by Estonian research 
projects SF0180013s07, IUT2-27 and by European Union through the Euro¬ 
pean Regional Development Fund (project 3.2.01.11-0029). 

A The effect of self-compression on the shape of transverse ILM 


Let us consider the effect of longitudinal displacements of atoms, caused by 
the transverse ILM, to its shape and frequency. From Eq. ([HD it follows that 
the TILM causes the change (reduction) of the atomic distances by Xc = 
—z/3 |/^(0)/4 i/i and slight vibrations of atoms in x direction with the frequency 
2uq (the forced vibrations). Neglecting the small vibrations mentioned, we hnd 
that on the average the x-dependent terms in Eq. (]HD give the following mean 
(averaged over period of transverse vibrations) contribution to the potential 
energy: 


t7 2 

= yx. 


+ = 


T28i/i 64i/i 


y\o) = 




1281/1 


y\0) 


(A.l) 


Let us hnd the nonlinear interaction Vc of a transverse motion with the mean 
value given by Eq. flA.ip . We restrict ourselves to the nonlinear interaction 
of oc type. Taking into account the relation = 3?/^(0)/8, we get W = 
—z/|?/^/48z/i. We suppose that this potential effectively accounts for the effect 
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of displacements of atoms in x direction. Adding this term to V(0,y) in Eq. 
(jS]), we get the following effective potential energy of the transverse vibrations 


Veff(xc, y) = V (0, y) + Vc = ^sy^ + ^ 2 /^, 


(A.2) 


where 1/4 = 1/4 — 2i>^/3ui. The TILM corresponding to this potential energy is 
described by Eq. ([6]) with the parameter 



(A.3) 


For the TILM with the amplitude Aq = 0.03 of the central atom in 5% 
stretched chain ^ 0.167. This gives e ~ 0.113 and Uq ~ 0.1883 in agreement 
with numerically found values of e and coq. 
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